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Abstract 


The spatial development of laminar-turbulent 
transition in a 45° swept-wing boundary layer is 
investigated using the large-eddy simulation ap- 
proach. Both stationary and traveling disturbances 
are initiated by steady and random forcings. The 
numerical simulations reproduce the surface shear 
streaks typical of the crossflow instability observed 
in experiments. Downstream of the transition loca- 
tion the wall shear stress is found to turn towards 
the streak direction. The strength of the stationary 
vortices grows exponentially, initially independent 
of the unsteady forcing amplitude and saturates at 
different levels during the later development, de- 
pending on the magnitude of the traveling vortices. 
Streamwise velocity contours show the evolution 
from a wave-like structure to the “half-mushroom” 
structure for the lower-amplitude, traveling-wave 
case, leading to double inflection points in the wall- 
normal velocity profiles prior to transition; the non- 
linear interactions involving the stationary vortices 
are also much stronger, as indicated by the ampli- 
fication of higher harmonics of the primary mode. 
The frequency spectrum of the traveling waves show 
a high frequency secondary instability prior to tran- 
sition that, however, is not prominent in the higher- 
amplitude traveling-wave case. 
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Introduction 

The transition process in a swept-wing bound- 
ary layer is significantly different and more complex 
than transition on a flat plate. One important type 
of instability that exists in the swept-wing bound- 
ary layer is the so-called crossflow instability, which 
is caused by the inflection point present in the sec- 
ondary (crossflow) velocity profile within the bound- 
ary layer. This instability plays a dominant role 
in the favorable pressure gradient region where the 
Tollmien-Schlichting waves are suppressed. 

The most distinctive feature of crossflow insta- 
bility is the existence of co-rotating stationary and 
traveling vortices, also known as crossflow vortices. 
The stationary vortices, which are initiated by sur- 
face roughness, tend to be aligned approximately in 
the streamline direction, and were observed in flow- 
visualization experiments (Dagenhart, 1992). Prior 
to transition, the velocity field is so strongly influ- 
enced by the stationary vortices that the stream- 
wise velocity profiles may exhibit multiple inflec- 
tion points (Dagenhart, 1992; Joslin and Streett, 
1994; Joslin, 1995). These inflectional profiles gen- 
erate higher-frequency, secondary-instability travel- 
ing modes (Kohama et al. , 1992). Experiments also 
showed that traveling vortices, especially the sec- 
ondary instability modes, play an important role in 
driving the laminar flow to transition and turbu- 
lence (Muller and Bippes, 1988). 

Despite many theoretical and experimental stud- 
ies, the current understanding of transition in a 
swept-wing boundary layer is still incomplete be- 
cause of the complexity of the problem, and because 
of the difficulties in carefully controlling the exper- 
imental conditions and the inclusion of nonlinear 
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effects in theoretical studies. Numerical simulations 
employing the full Navier- Stokes equations — di- 
rect numerical simulations (DNS) — can, therefore, 
give a substantial contribution to the understand- 
ing of the flow physics because the ingestion of dis- 
turbances into the computed boundary-layer can be 
more carefully controlled than in experiments. The 
cost of DNS limits, however, its application to flows 
of engineering interest at high Reynolds number. 
Even at moderate or low Reynolds numbers, the 
computational effort required to perform a calcula- 
tion of the full transition from laminar to turbulent 
flow is extremely high. 

Large-eddy simulation (LES) is an alternative, 
less expensive, approach for the numerical solution 
of turbulent and transitional flow problems. Un- 
like DNS in which all scales of motion are resolved, 
only the dynamically important large scales need 
to be resolved in LES, which results in a signifi- 
cant decrease of the CPU time required for a sim- 
ulation compared with DNS. The influence of the 
unresolved small scales on the large scales is mod- 
eled. Since the small scales tend to be more homoge- 
neous and isotropic than the large scales (which are 
highly flow- dependent), the modeling can be simpler 
and more universal than that used in the Reynolds- 
averaged approach. 

The objective of this work is to study of transi- 
tion in a swept- wing boundary layer using the large- 
eddy simulation approach. Two simulations are 
conducted with small and large initial amplitudes 
of the time-dependent disturbances and fixed am- 
plitudes of the stationary disturbances. The effect 
of random, three-dimensional, and time-dependent 
disturbances on the flow breakdown, and the in- 
teraction between the stationary vortices typical of 
swept-wing transition and the traveling modes gen- 
erated by the time-dependent disturbance will be 
investigated to elucidate the mechanisms by which 
laminar flow becomes turbulent. 

Numerical Method 

In the present study, the filtered governing equa- 
tions for the disturbance variables (iq, p) are in the 
form: 


and 


where an overbar denotes the application of the 
filtering operation, and Tij = lljUj — uplj is the 
subgrid-scale (SGS) stress tensor, which must be 
modeled. The base flow field (U 0 i) is given by 
the Falkner-Skan- Cooke similarity solution (Cooke, 
1950). The Reynolds number is Re = Q r 8/v , where 
Q r is the local edge velocity at the computational 
inflow boundary, and S = ^/ra c /U oe is the reference 
length at the inflow. 

The filtered Navier- Stokes equations are solved 
using the fractional time step method (Chorin, 
1968). Fourth-order finite difference and fourth- 
order compact difference schemes are used in the 
streamwise direction (x) for the pressure and the 
momentum equations, respectively, Chebychev se- 
ries are used in the wall-normal direction (y), and 
Fourier series in the spanwise direction (z). The 
implicit Crank-Nicolson time- advancement is used 
for the wall-normal diffusion terms and a three- 
stage Runge-Kutta scheme for the remaining terms 
(Joslin et al , 1992). The subgrid-scale stress term 
is modeled by a localized version (Piomelli and Liu, 
1995) of the dynamic SGS model (Germano et al , 
1991). Spanwise averaging as well as three point 
averaging in the streamwise and wall-normal direc- 
tions are used when calculating the localized dy- 
namic SGS model coefficient. The total viscosity 
(molecular + eddy viscosity) is also forced to be 
non-negative to ensure numerical stability. 

One of the major difficulties associated with 
the numerical simulation of spatially developing 
boundary-layer transition is to define the outflow 
boundary conditions. Here the buffer domain tech- 
nique proposed by Streett and Macaraeg (1989) is 
used, in which the governing equations are gradually 
parabolized in a buffer region which is appended at 
the end of the computational domain, thus eliminat- 
ing the necessity of applying the outflow boundary 
conditions. A fourth-order compact-difference fil- 
ter (Lele, 1992) is used in the streamwise direction 
every 10 time steps, to eliminate numerical wave 
reflections at the outflow (Pruett et al . , 1995). 

Simulation Parameters 
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In the present simulations, some of the parame- 
ters were chosen to match the DFVLR 45° swept 
flat-plate experiment carried out by Muller and 
Bippes (1988). The base flow was determined by 
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imposing the sweep angle and the pressure distribu- 
tion 

c p (x c ) = 0.941 — 0.845x c (3) 

on the Falkner-Skan- Cooke equations. Two compu- 
tational boxes with different grid resolutions were 
used to cover the whole physical domain to be stud- 
ied, and the inflow conditions for the second box 
were read from the outflow of the first one. This 
implementation was used successfully in the large- 
eddy simulation of subharmonic breakdown in a flat- 
plate boundary layer (Huai, 1996; Huai et al , 1997). 
Near the computational inflow, small disturbances 
were introduced by imposing a wall-normal velocity 
distribution in the form of a suction and blowing 
strip. A dominant stationary vortex mode was ex- 
plicitly introduced by applying steady suction and 
blowing. The dominant spanwise wavelength was 
found to be = 36 for the conditions used in 
Joslin and Streett (1994), and the same conditions 
and dominant wavelength were used in the present 
study. Traveling crossflow vortices were introduced 
into the flow by random modulations of the suction 
and blowing at each time step. The zero- frequency 
mode was filtered out from the random part of the 
suction and blowing to prevent the traveling-wave 
generating mechanism from altering the initial sta- 
tionary wave, and also to maintain overall mass con- 
servation. Since no other restriction was applied, 
the random function contained all the non-zero fre- 
quency and wavenumber modes allowed by the grid 
resolution, time step size and time interval of the 
data sampling. 

Two simulations, denoted by LOW-T and HIGH- 
T, were carried out. While the steady suction and 
blowing had the same magnitude in both simula- 
tions (A 0 = 10 — 4 , where A 0 is the amplitude of 
the wall-normal velocity fluctuation in the suction 
strip), the magnitudes of the unsteady suction and 
blowing, At, were different, with A t = 10 -4 and 
1(T 2 for the LOW-T the HIGH-T cases, respec- 
tively, to study the role of traveling vortices and 
their effects on the transition process. The Reynolds 
number at the inflow was 998, which is equivalent 
to a chord Reynolds number Re c = 2.73 x 10 6 . The 
two simulations covered the chordwise regions of 
0.2 < x c < 0.54 and 0.2 < x c < 0.5, and required 
about 100 and 90 Cray C-90 hours respectively. 
The two boxes for LOW-T had streamwise x wall- 
normal x spanwise lengths of 358x50x108 and 
386 x 50 x 108, and used 326 x 49 x 33 and 461 x 61 x 49 
grid points. The two boxes for HIGH-T had do- 


main lengths of 272 x 50 x 108 and 386 x 50 x 108, and 
used 266x49x33 and 461x61x49 grid points re- 
spectively. 

Results 

Code validation 

The numerical scheme and SGS models had previ- 
ously been validated in two studies: the first was the 
calculation of subharmonic instability in a flat-plate 
boundary layer (Huai, 1996; Huai et al , 1997); the 
second, the calculation of transition in a swept wing 
in presence of stationary disturbances only (Huai et 
al , 1994). 

The first calculation showed good agreement be- 
tween the LES results and linear theory in the 
early transition stages, and excellent agreement 
with the data obtained from DNS calculations 
(Spalart, 1986). The second simulation studied a 
flow very similar to the one under consideration 
here; the results were compared with the DNS data 
of Joslin and Streett (1994), and good agreement 
was achieved between the LES and DNS data. 

A thorough grid-resolution study could not be 
performed due to the computational expense in- 
volved; however, one test was run using the same 
conditions as case HIGH-T, but a refined mesh in 
the second box (three times as many points in x, and 
twice as many in z; the spanwise domain was one- 
third of that used in simulation HIGH-T, however). 
The fine-mesh calculation was in good qualitative 
agreement with the coarser one; for x c < 0.42, the 
agreement was also quantitatively satisfactory. 

A direct comparison of the present LES results 
with experimental or theoretical data was not per- 
formed in this study, given the wide variety of flow 
conditions used in the previous studies. Several flow 
features obtained from the LES agree well with ex- 
perimental and theoretical data (see below). Based 
on the detailed comparisons of the LES results with 
DNS for flat-plate boundary-layer transition and the 
grid refinement comparisons for the present three- 
dimensional flow reported by Huai (1996), it is pos- 
sible to conclude that the early nonlinear stages 
of transition are captured quite accurately, while 
the marginal grid resolution in the final breakdown 
stages allows us to draw only qualitative conclusions 
regarding the transition itself. 

Flow visualization 

In flow visualization experiments of three- 
dimensional boundary-layer transition, surface 
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Figure 1. Contours of mean wall shear 
log 10 (dV/dy) w . Top: LOW-T; bottom: HIGH-T. 
The white vertical lines represent the junction be- 
tween the first and the second boxes. 


streaks, which are the signatures of stationary cross- 
flow vortices, have been frequently observed. Poll 
(1985) was able to observe the surface streaks in a 
swept cylinder boundary layer. Dagenhart (1992) 
documented the stationary crossflow vortices in a 
swept-wing boundary layer by taking photographs 
of the naphthalene flow visualization: a white naph- 
thalene sublimating coating was placed over the 
black swept-wing surface; the existence of station- 
ary crossflow vortices gave rise to alternating high- 
and low-shear regions on the wall. 

Contours of mean wall shear from the present sim- 
ulations are shown in Fig. 1 (two computational do- 
mains are shown side by side in each figure). The 
alternating high- and low-shear regions are evident. 
As in the experiment of Dagenhart (1992), no vor- 
tex merging or adjustment of vortex spacing is ob- 
served, in contrast to the experimental results of 
Arnal and Juillen (1987). The vortex adjustment 
in the numerical simulations could be restricted to 
some extent by the spanwise periodicity. 

The transition location can be estimated in Fig. 1 
as the place where the distinct streak patterns be- 
come blurred or disappear. For LOW-T and HIGH- 
T they are at x c ~ 0.45 and 0.42 respectively. Tran- 
sition occurs earlier in the higher time-dependent 
disturbance (HIGH-T) case. This result is consis- 
tent with the experimental results of Muller and 
Bippes (1988), who suggested that the traveling 
waves play the most important role in determining 
the onset of transition, while large amplitudes of the 
stationary vortices do not seem to advance transi- 
tion. Another interesting observation is that in the 


LOW-T case the change of the streak patterns in 
the downstream region is milder and remnants of 
the stationary vortices extend deep into the transi- 
tion and turbulent regions. 

Development of the mean quantities 

Time- and spanwise- averages were computed 
from the sampled flow-field data to obtain mean- 
flow quantities. Instead of using the velocity 
components in the body-oriented coordinate frame 
(x,y,z), the velocities were projected into the flow- 
oriented coordinate system ( x s ,y s ,z s ) defined by 
the local external streamwise (x s ) and the cross- 
flow (z s ) directions. The displacement thickness 5*, 
momentum thickness 0 S , shape factor H s and skin- 
friction coefficient C/ jS associated with the stream- 
wise mean velocity U s are defined in this coordinate 
system as 


* - L <4) 



where U s ^ e is the boundary-layer edge velocity. 

Figure 2 shows the chord wise variation of the 
skin-friction coefficient and shape factor for the 
LOW-T and HIGH-T cases. Initially, while the 
disturbances are growing, the mean flow is unal- 
tered from its laminar state. As the crossflow dis- 
turbances amplify, however, the mean flow begins 
to be affected by the disturbances, the skin fric- 
tion increases and the shape factor decreases from 
their laminar values, indicating the onset of transi- 
tion to turbulence in imminent. Finally the curves 
level off to a turbulent state. The higher time- 
dependent perturbation case, HIGH-T, results in 
earlier transition than LOW-T, but the differences 
between the shape factors are small in the early 
transition regions and become pronounced only in 
the later stages of transition. The laminar shape 
factor is lower than the flat-plate boundary-layer 
value of 2.6, because of the three-dimensionality of 
the boundary-layer flow. The turbulent values are 
almost the same for both LOW-T and HIGH-T, re- 
gardless of the difference in A t . 
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Figure 2. Chordwise variation of the (a): the skin 
friction coefficient C/ jS , and (b): shape factor H s . 
: LOW-T; : HIGH-T. 

Due to both the sweep angle of the flow and the 
pressure gradient, the streamlines and wall shear- 
stress change direction with chordwise location. As 
shown by Joslin and Streett (1994), the crossflow 
vortices follow the external streamwise direction, 
within a few degrees. The angles obtained in the 
present study by measuring the streak patterns 
are plotted in Fig. 3. The deviation, e — 4°, 
from the direction of the external streamwise di- 
rection f3 matches the experimental result of e — 
3° ~ 5° (Nitschke-Kowsky and Bippes, 1988). The 
wall shear-stress initially follows the laminar base- 
state direction; as transition occurs, however, the 
strength of the crossflow disturbances affects the 
near- wall flow and the shear direction switches from 
the laminar state to the direction of the crossflow 
vortex disturbance. 

Development of the stationary vortices 

In Fig. 4, the maximum streamwise velocity as- 
sociated with the stationary vortices for both sim- 
ulations is shown. Irrespective of the different 
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Figure 3. Chordwise variation of velocity, vortex 
and wall-shear angles. External streamwise angle 
(3: ; stationary crossflow vortex angle f3 — e, 


LOW-T: □ , HIGH-T: A ; f3 - 4°: ; wall shear- 

stress angle of the laminar base flow: ; actual 
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Figure 4. Chordwise growth of the stationary dis- 
turbances. LOW-T; HIGH-T. 


amplitudes of the traveling vortices, both simula- 
tions give stationary disturbances with exponential 
growth rates oti — 11. This growth rate is very close 
to the experimental result, oti — 13 ± 1.3, of Muller 
and Bippes (1988) for the Tu— 0.05% case, given the 
experimental uncertainty. The experimental growth 
was also in good agreement with results (a* = 12) 
of a numerical calculation (Meyer, 1987). However, 
the experiment also showed that the initial ampli- 
fication of the stationary vortices was lower for the 
higher freestream turbulence, i.e., Tu— 0.15%, case, 
apparently in contrast with the present results, that 
yield the same growth rate for both HIGH-T and 
LOW-T cases. In the Tu— 0.15% case, however, 
the initial magnitude of the traveling disturbance 
was about the same as that of the stationary one, 
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and about half that of the stationary disturbance 
for the low turbulence case. In the present simula- 
tions, on the other hand, the traveling disturbance 
is weaker, its amplitudes is close to the experimen- 
tal Tu = 0.05% case in HIGH-T, and one order 
of magnitude smaller than the stationary wave; in 
LOW-T, the traveling wave is three orders of magni- 
tude smaller than the stationary wave. The present 
results, combined with the experimental findings, 
seem to suggest that if the magnitude of the trav- 
eling waves has not exceeded a certain threshold 
(somewhere between 0.5 and 1 of the amplitude of 
the stationary waves), the initial amplification of the 
stationary crossflow vortices can be predicted satis- 
factorily by the theory (Meyer, 1987); otherwise it 
will be lower than the theoretical prediction. The 
reason for this is that the accelerated nonlinear in- 
teractions transfer energy from the stationary waves 
to the traveling ones. 

The saturation (or over-saturation for LOW-T) of 
the stationary vortices prior to transition can also 
be observed in Fig. 4. It starts at different chordwise 
distances depending on the locations of transition. 
This observation is in agreement with the findings 
of the crossflow instability experiment of Muller and 
Bippes (1988), which showed that earlier transition 
in the case of stronger traveling waves results in ear- 
lier saturation. 

Figure 4 shows the growth of mean distortion of 
the laminar base flow, but it provides no informa- 
tion on the wavenumber content of the stationary 
vortices, their growth, and on the possible non- 
linear wave interactions that may play an impor- 
tant role in understanding the transition process. 
For example, in the crossflow experiment of Saric 
and Yeates (1985), the stationary vortices having 
the wavelength of 1 cm predicted by linear theory 
were observed, while in the hot-wire measurements 
the second harmonic having a wavelength of 0.5 cm 
away from the wall was also found, and its charac- 
teristics amplitude was as high as three times that of 
the 1 cm primary wave. Reed (1987) studied wave 
interactions in the framework of secondary instabil- 
ity and concluded that the anomalies found in the 
experiment are introduced by a crossflow/crossflow 
interaction between the primary vortices and non- 
stationary disturbances whose wave-length is half 
that of the primary crossflow vortices. 

To investigate the nonlinear interactions that oc- 
cur in this flow, the spanwise wavenumber spectra 
of the steady streamwise perturbation velocity were 
calculated. Figures 5 and 6 give the chordwise vari- 



Figure 5. Chordwise development of the spanwise 
wavenumber spectrum for the LOW-T case. 



Figure 6. Chordwise development of the spanwise 
wavenumber spectrum for the HIGH-T case. 
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ation of maximum magnitude for each mode defined 
as 

Amax(%->kz) = max [A(x, y, A^)], 

V ^ [Oj Vmax\ (8) 

where A(x,y,k z ) is the magnitude of mode k z at 
point (x, y). The wavenumber of the primary cross- 
flow mode (denoted by K z ) satisfies the relation- 
ship K z L z /2tt = 3, since the steady suction/blowing 
has three periods over the spanwise box length. 
In the lower-amplitude traveling-wave case the pri- 
mary crossflow mode reaches an amplitude equal to 
28% of the freestream velocity prior to transition, 
while the second harmonic 2 K z attains the peak 
magnitude of 18%. It is also evident in the figure 
that the appearance of higher harmonics lags behind 
that of lower ones, and that their magnitudes de- 
crease as the wavenumber increases. This suggests 
that the growth of higher harmonics is not due to 
the secondary instability analyzed by Reed (1987), 
but to the nonlinear wave interactions between the 
stationary vortices themselves (the second harmonic 
is excited by the interaction between the primary 
wave and itself, the third harmonic by the the pri- 
mary mode and the second harmonic, etc.) The 
saturation of the stationary vortices must also be 
caused by these nonlinear wave interactions, since 
the amplitude of the traveling waves is still small 
at this stage. In fact, Huai (1996) found a higher 
growth rate of secondary instability for the traveling 
modes, since the stationary vortices have reached fi- 
nite amplitudes. In the higher- amplitude traveling- 
mode case (Fig. 6), however, the higher harmon- 
ics never amplify to the same magnitude, and the 
primary wave itself only reaches about 14% of the 
freestream velocity, due to early nonlinear wave in- 
teractions of stationary /traveling type. 

The chordwise growth of several harmonics and of 
the spanwise-averaged perturbation velocity, the ze- 
roth wavenumber mode, are shown in Fig. 7. In the 
LOW-T case (Fig. 7a) the fundamental mode grows 
exponentially (a* = 11) and is the largest one be- 
cause the steady suction and blowing generates only 
this single mode; the second harmonic, 2 K z , has a 
higher amplification rate (a* = 24), the third har- 
monic, 3 K z , generated by the fundamental interact- 
ing with the first harmonic mode, has an even higher 
amplification rate (a* = 39). At x c = 0.39, the 
fundamental mode starts to saturate and reaches a 
peak amplitude 28% at x c = 0.41; the second har- 
monic saturates at approximately the same location, 




Figure 7. Chordwise development of several station- 
ary modes, : mean velocity; : K z mode; 

: 2 K z mode; : 3 K z mode, (a): LOW-T 

case; (b): HIGH-T case. 


while the saturation of the third harmonic occurs 
downstream of the first- and second-mode satura- 
tion point. The mean velocity mode also shows the 
tendency to saturation when the primary mode sat- 
urates, but it regains exponential growth when the 
transition zone is entered. At the transition loca- 
tion x c = 0.45, all the modes have amplitudes of 
the same order of magnitude. In the HIGH-T case 
(Fig. 7b) the primary mode and the second har- 
monic have the same initial growth rates as in the 
LOW-T case, and begin to saturate at x c = 0.36, 
while the exponential growth of the third harmonic 
is sustained further downstream to x c = 0.39. The 
higher harmonics are one order of magnitude smaller 
than the primary mode at saturation, and they are 
not evident in Fig. 6. The growth of the zeroth 
wavenumber follows that of the second harmonic 
very well up to the saturation in both cases. 

Figure 8 shows the streamwise velocity contours 
in the yz-pl&nes parallel to the leading edge of the 
swept wing. At x c = 0.3, the crossflow vortices are 
still very weak and the mean flow is uniform in the 
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Figure 8. Contours of the time-averaged streamwise 
velocity U s for the LOW-T case, (a) x c = 0.3; (b) 
x c = 0.35; (c) x c = 0.4; (d) x c = 0.45; (e) x c = 0.5; 
(f) Xc = 0.54. 


spanwise direction. The presence of the crossflow 
vortices becomes detectable near x c = 0.35, and the 
low and high wall-shear regions become very clear at 
x c = 0.40. When the transition location x c = 0.45 
is reached (Fig. 8d), the distortion of the mean ve- 
locity is so strong that the so-called “half mush- 
room” structures are observed (Dagenhart, 1992; 
Joslin and Streett, 1994). This structure (or vortex 
roll-over state) is accompanied by double inflection 
points in the wall-normal profiles of the streamwise 
velocity. Besides the primary alternating wall-shear 
pattern, weak secondary crossflow vortices can also 
be observed relatively close to the wall at z & 10, 45 
and 80. Malik et al. (1994) also observed these sec- 
ondary vortices in their numerical study. As the flow 
goes into the transition and turbulent regions, the 
strength of the vortical structures gradually dimin- 
ishes, but the spanwise periodicity remains observ- 
able; this explains why the streaks extend into the 
downstream region in the wall-shear plot (Fig. 1). 

For the higher-amplitude traveling-wave case 
(HIGH-T) the streamwise velocity contours (Fig. 9) 
show the spanwise distortion of the mean flow, but 
the “half mushroom” structure that was observed 



Figure 9. Contours of the time-averaged streamwise 
velocity U s for the HIGH-T case, (a) x c = 0.3 (b) 
x c = 0.35; (c) x c = 0.4; (d) x c = 0.42; (e) x c = 0.45; 
(f) x c = 0.5. 

in the LOW-T case is not present because the sta- 
tionary vortices are not as strong. Correspondingly, 
the wall-normal velocity profiles in the HIGH-T case 
exhibit only marginal inflection points. 

Frequency analysis 

It has been shown by linear stability theory and 
experiments (Mack, 1984; Muller and Bippes, 1988; 
Dagenhart, 1992) that the crossflow velocity pro- 
files present in three-dimensional swept-wing flows 
can excite both stationary and traveling vortex 
disturbances, and that the traveling waves have 
higher growth rates than the stationary ones. In 
low-turbulence environments such as those encoun- 
tered in flight conditions, however, receptivity the- 
ory (Crouch, 1993) indicates that the stationary dis- 
turbances will have much larger amplitudes com- 
pared with the traveling modes because of the low 
free-stream turbulence levels. However, Muller and 
Bippes (1988) note that the traveling waves play 
a more important role than the stationary waves 
in determining the onset of transition. Although 
the frequency range of the most amplified primary 
traveling waves can be predicted satisfactorily by 
the linear stability theory, Kohama et al (1991) 
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Figure 10. Frequency spectrum of the streamwise 
velocity u s at the maximum traveling wave locations 

for the LOW-T case. x c = 0.3, x c — 0.35, 

x c — 0.4, x c — 0.45, x c — 0.5, 

x c = 0.54, + f~ 5 / 3 . 

were able to detect a most amplified high-frequency 
band prior to transition in a swept wing boundary- 
layer experiment in the ASU low-turbulence wind 
tunnel (Tu= 0.02%) that could not be predicted by 
the linear stability theory. This was believed to 
be due to the secondary instability associated with 
the strongly distorted streamwise velocity profiles. 
The same phenomena was reported earlier with- 
out explanation by Poll (1985) for a swept-cylinder 
boundary-layer experiment. 

This issue is further studied by examining the de- 
velopment of the frequency spectra of the stream- 
wise fluctuation velocity u s . For the LOW-T case 
(Fig. 10), a low frequency band centered at / = 
4.17xl0 -3 begins to amplify, in agreement with the 
prediction of linear stability theory. At x c = 0.35, 
however, a second peak is observed in the spectrum, 
corresponding to the amplification of the second 
harmonic of the primary traveling wave. The travel- 
ing disturbance is dominated by the low-frequency 
band up until x c = 0.4, but at x c = 0.45 a high- 
frequency band with frequencies one order of magni- 
tude higher is strongly amplified, in agreement with 



Figure 11. Frequency spectrum of the streamwise 
velocity u s at the maximum traveling wave locations 

for the HIGH-T case. x c = 0.3, x c = 

0.35, x c — 0.4, x c = 0.42, x c = 

0.45, x c = 0.5, + /~ 5 / 3 . 

the findings of Kohama et al. (1991), and the peak in 
the spectrum moves to the second harmonic. In the 
transition and turbulent regions, the distribution of 
the spectrum becomes typical of a turbulent flow, 
exhibiting the —5/3 slope of the inertial subrange. 
It should be pointed out that due to the grid resolu- 
tion and accuracy of the numerical scheme, a drop 
of about four orders of magnitude can be obtained 
in the spectrum; very small values of the spectrum 
are produced by either aliasing error or power leak- 
age from higher energy modes and are not physical. 
However, the LES modeling of the small-scales is 
producing credible spectral trends. 

Figure 11 shows the frequency spectrum for the 
HIGH-T case. A shift of the peak low-frequency 
mode occurs at x c = 0.4. However, the high fre- 
quency band is not prominent at the transition lo- 
cation (x c = 0.42); the secondary instability asso- 
ciated with the inflectional points is much weaker 
because the stationary disturbance is weaker and 
exhibits only marginal inflection points. Because 
of the earlier nonlinear stationary /traveling wave 
interaction the spectrum contents are higher and 
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much closer to those of turbulent flows. 

Conclusions 

The spatial development of crossflow instability in 
a 45° swept-wing boundary layer was studied using 
the large-eddy simulation (LES) technique. The nu- 
merical simulations successfully reproduce the wall- 
shear streaks typical of the crossflow disturbance. 
The simulations also show that transition occurs 
later in the simulation with weaker traveling waves, 
as expected, and that the remnants of the streaks 
extend deeply into the transition and turbulent re- 
gions in the low unsteady-forcing case. The tran- 
sition locations obtained from the streak patterns 
compare well with those obtained from the varia- 
tion of the streamwise skin-friction coefficient. The 
streaks are approximately aligned in the direction 
of the external streamline, forming a 4° differential 
angle with the freestream, in good agreement with 
experimental results. Downstream of the transition 
location the wall shear stress is found to align with 
the streak direction. 

The strength of the stationary vortices is found 
to grow exponentially, initially independent of the 
unsteady forcing amplitude; however, it saturates 
at different levels during the later development, de- 
pending on the magnitudes of the traveling vor- 
tices. Streamwise velocity contours in the yz-pl&nes 
show the evolution from a wave-like structure to the 
so-called “half mushroom” structure for the lower- 
amplitude traveling wave (LOW-T) case, leading to 
double inflection points in the wall-normal veloc- 
ity profiles prior to transition. However, the “half- 
mushroom” structure is not observed in the HIGH- 
T case. The nonlinear interactions involving the 
stationary vortices are much stronger in the low 
unsteady-forcing case, as indicated by the ampli- 
fication of higher harmonics of the primary mode. 

The numerical results also show that the traveling 
waves are dominated early by low frequency com- 
ponents. Prior to transition, in the low unsteady- 
forcing case a secondary instability band with high 
frequencies is most amplified, followed by rapid 
spectrum filling. However, this secondary instabil- 
ity is much weaker in the higher unsteady-forcing 
case. 
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